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Background: Species in the Culex pipiens complex are common almost all over the world and represent 
important vectors for many serious zoonotic diseases. Even if, at the moment, many of the pathogens 
potentially transmitted by Cx. pipiens are not a problem in northern Europe, they may, with increasing 
temperatures and changing ecosystems caused by climate change, move northward in the future. Therefore, 
the question whether or not the Cx. pipiens populations in northern Europe will be competent vectors for 
them is of high importance. One way to estimate the similarity and the rate of contact between European Cx. 
pipiens populations is to look at the gene exchange between these populations. 

Methods: To test the genetic diversity and degree of differentiation between European Cx. pipiens 
populations, we used eight microsatellite markers in 10 mosquito populations originating from northern, 
central, and southern Europe. 

Results: We found that three of the analyzed populations were very different from the rest of the populations 
and they also greatly differed from each other. When these three populations were removed, the 
variance among the rest of the populations was low, suggesting an extensive historic gene flow between 
many European Cx. pipiens populations. 

Conclusions: This suggests that infectious diseases spread by this species may not be associated with a certain 
vector genotype but rather with suitable environmental conditions. Consequently, we would expect these 
pathogens to disperse northward with favorable climatic parameters. 
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Mosquitoes are found throughout the world 
except in locations that are permanently frozen 
(1, 2). In almost all mosquitoes, the female 
feed on vertebrate blood to obtain the protein they need 
for egg development. A complex salivary secretion in 
mosquitoes helps the feeding process, and it is the 
injection of this fluid to the host capillaries that enables 
several kinds of simple organisms - viruses, protozoa and 
nematodes - to enter the vertebrate blood system (1, 2). 
This process makes it possible for pathogens to 
exploit mosquitoes as a principal vector of infectious 
diseases in both developing countries and in industrial 
nations (3). In many, mainly tropical and economically 
challenged parts of the world, mosquito-borne diseases 



are flourishing and they are responsible for a large 
proportion of the global morbidity and mortality (2). 

Species in the Culex pipiens complex are common 
almost all over the world and act as important vectors for 
many serious zoonotic diseases, mostly because they 
often are the most common mosquitoes in urban areas 
(4). The taxonomy and systematics of the Cx. pipiens 
complex is a challenge. The form molestus that was 
originally described as a species by the Linnaean disciple 
Forsskal in 1775 from Egypt has since been reduced to 
represent a bioform (bf.) of the species Cx. pipiens. This 
decision was based upon the fact that no morphological 
characters were found that would allow the discrimina- 
tion of one species from the other. Fonseca et al. (4) 
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showed that there were distinct genetic differences using 
microsatellites between populations of bf. pipiens and 
molestus and that they probably had undergone a recent 
speciation event. Because bf. pipiens is strictly ornitho- 
philic while molestus readily bites humans and other 
mammals, this is of epidemiological importance when 
discussing diseases such as West Nile fever (WNF). 
Supposedly, molestus will act as a bridge vector between 
the bird reservoirs and susceptible mammal hosts, while 
pipiens only will maintain the bird-to-bird transmission. 
Except for some underground populations of Cx. pipiens 
bf. molestus, in northern European cities, most specimens 
of this temperate species require a blood meal for egg 
development (4). 

Anthropogenic climate change is presently exerting 
multifaced effects on ecosystems, communities, and 
populations (5). The ecology, behavior, development, 
and survival of mosquitoes and the disease transmission 
dynamics are strongly influenced by climatic factors. 
Temperature, rainfall, and humidity are thought to be 
the main components influencing the life history of 
mosquitoes, even if other factors, such as the wind and 
photoperiod can have significant influence (1, 6). The 
same factors also influence the survival and transmission 
pattern of mosquito-borne pathogens - for example, the 
temperature affects the rate of pathogen multiplication 
within the insect, which in turn affects the rate of infective 
organisms in the salivary secretion and thus the likelihood 
of successful transmission to another host (1, 2, 7). 

Understanding the connection between infectious dis- 
eases and climate is difficult because of the multivariate 
nature of climate change and non-linear thresholds in 
both disease and climate change processes. It is commonly 
assumed that the distribution of infectious diseases is 
determined by climate and that the present climate change 
will increase their incidence and geographic distribution. 
Consequently, the relationship between temperature and 
the transmission of pathogens has gained substantial 
attention (8). However, just because a vector species can 
transmit a certain pathogen in one geographical area, that 
does not necessarily mean that the populations of the 
same species are competent hosts for the same pathogen in 
another geographic location (9) as the vector competence 
of mosquitoes can be influenced by both temperature and 
viral genetics (8). 

A wide variety of pathogens can be transmitted by Cx. 
pipiens. For example, these mosquitoes are important 
vectors of West Nile virus that was first characterized in 
Uganda in 1937 but has spread over large areas of Africa, 
Europe, the Middle East, Central Asia, and the Americas. 
Usutu virus is another example of a previously tropical 
mosquito-borne virus that has now been introduced and 
established in temperate climate - originally isolated from 
a bird-biting mosquito, Culex univittatus, in South- 
Africa, and historically detected in various places in 



sub-Saharan Africa, the virus, in 2001, suddenly emerged 
in the Austrian capital of Vienna (10) and recently in Italy 
(11, 12) and in southern Germany (13). Dirofilariases are 
important vector-borne parasitic zoonoses that mainly 
affect the canine population. Cx. pipiens has been sug- 
gested to be one of the main vectors for the heartworm 
disease by Dirofilaria immitis that can be rather common 
in some areas in southern Europe (14). Rift Valley fever 
virus is another virus that can be carried by Cx. pipiens, 
and can cause severe disease in both animals and humans 
in various parts of Africa (however, carried by other 
mosquito species in African cases) (15). In September 
2000, Rift Valley fever cases were diagnosed in Saudi 
Arabia and Yemen, raising concerns that it could extend 
to other parts of Asia and Europe. 

The above-mentioned pathogens have not yet been 
causing problems in northern Europe, however with 
increasing temperatures and changing ecosystems caused 
by climate change, it is possible that they continue to 
move northward. Therefore, the question whether or not 
the Cx. pipiens populations of northern Europe will be as 
competent vectors for mentioned pathogens as the more 
southern populations is of high importance for risk 
analysis. 

One way to estimate the similarity and the rate of 
contact between European Cx. pipiens populations is to 
look at the gene exchange between mosquito populations. 
Since gene flow is always lower than the actual movements 
of individuals, this will give a rather conservative estimate 
of the amount of spread of individuals between parts of 
Europe. By analyzing the extent of contact between 
populations that have been shown to be able to carry 
certain pathogens and populations in northern Europe 
that have yet not been seen to transmit them, we can 
estimate the possibility of natural pathogen dispersal 
under changing climatic conditions. To estimate the 
genetic diversity and degree of differentiation between 
European Cx. pipiens populations, we used eight micro- 
satellite markers in 10 mosquito populations originating 
from northern, central, and southern Europe. 

Methods 
Sampling 

Summer 2009, adult individuals of Culex pipiens were 
caught by using different models of C0 2 and light traps in 
rural habitats at 10 different European locations (Fig. 1). 
Mosquitoes were preserved in 90% ethanol and sent to 
Uppsala by mail by the local mosquito researchers. Each 
Cx. pipiens population was represented by 30 adult 
individuals, all from wild, over-ground populations. 

Molecular analysis 

DNA from individual mosquitoes was extracted by a 
Chelex procedure proteinase-K-based extraction (16) 
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Fig. 1. Geographic origin of mosquito populations. 



using the protocol provided by the manufacturer. 
We used eight primers developed for Cx. pipiens (4, 17- 
19) (CQ13, CQ29, CxpGT9, CxpGT20, CxpGT40, 
CxpGT46, CxpGT51, and CxpGT53) and followed the 
PCR-protocols given in the original paper. The PCR- 
products were run on a MegaBase using dyed primers. 

Statistics 

We used the software Microchecker (20) to check for null 
alleles and other technical problems. We detected a 
significant number of null alleles in several loci, apparent 
as an uncommonly large proportion of homozygotes for a 
few alleles, which created large deviations from Hardy- 
Weinberg equilibrium. The repeatability of the amplifica- 
tions was 95%; so, there are no reasons to believe that 
technical problems caused this pattern. Since several 
samples did not amplify at all for a given locus, there 
are reasons to suspect that there are alleles present that 
did not amplify. Null alleles are commonly found in 
Lepidoptera (21) and seem to be present also in this 
Dipterian sample. The main effect of null alleles seems to 
be an overestimate of F ST (22), but no major distortions 
of genetic structure. We used the method by Chapuis and 
Estoup (23) to correct for the presence for null alleles 
since this method performed best in a comparative study. 
The method is implemented in the software FreeNA (23). 

To perform basic statistics, such as the number of 
alleles, allelic richness, and expected heterozygosity, we 
used the softwares Arlequin 3.5.1.2 (24) and Fstat 2.9.3.2 
(25). To calculate pairwise Fs T , we used FreeNA using 
both uncorrected data and data corrected for null alleles; 
95% confidence intervals were obtained through 1000 
bootstrap replicates. Since F ST is dependent on the level of 
heterozygosity, a low estimate can be due to a high within 



population variation rather than lack of differentiation. 
We therefore also calculated G' S t (26), which rescales 
F ST to have bounds between zero and one to get an 
intuitive picture of the degree of differentiation. To tests 
for possible isolation by distance patterns, we used the 
program Mantel 2.0 (27). Finally, to make neighbor- 
joining (NJ) trees based on the distances obtained, we 
used MEGA5 (28). To assess if the sampling and the 
number of loci used was enough to detect a structure if 
there is one (i.e. power of the tests), we used the program 
POWSIM (29). We used all loci and simulated i^-values 
for two populations being separated for 250 generations, 
each with an effective population size (N e ) of 5000 and 
with a sampling effort of 30 individuals (as in this study). 
The values were used to simulate a case where the split of 
the populations are very recent, and the populations are 
very large. Hence, this is a very conservative test of power, 
as if time is longer and N e is smaller, it would be much 
easier to find significant results. Nevertheless, even with 
these very conservative start values, we found that the 
expected F ST is 0.025 (95% range, 0.0139-0.0392) and 
the power 0.99. Thus, we are confident that, given the 
number of loci and sampling effort, we can find differ- 
ences if they exist. 

Because of some possible technical errors, we had to 
adjust our data. Consequently, we were not able to 
estimate genetic differentiation, as our data adjustments 
cannot cope with functions such as, for example, the 
estimation of gene flow. When the data are adjusted for 
null alleles (alleles that do not amplify), multilocus 
estimates of gene flow are no longer possible because 
the different loci are now made independent. Instead we 
used so-called private alleles (i.e. alleles only found in one 
population; 30) assuming that there were a number of 
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private alleles that did not amplify (null alleles). With help 
of GenePop and by using private alleles, we estimated the 
number of migrants per generation across all populations. 
The aim of estimating the gene flow with this method was 
to test if the observed amount of divergence was affected 
by gene flow or if it was what we could expect when 
no gene flow has occurred since separation. To implement 
this, we calculated expected i<sT by simulations using 
POWSIM with N c values ranging from 100 to 10,000 and 
time since separation (t) ranging from 100 to 5000 
generations and allele frequencies taken from the ob- 
served data. This will give us a range of expected F ST - 
values for different combinations of N c and t (for results, 
see Fig. 4). Populations from Romania, Turkey, and 
France were excluded from this calculation since they 
probably belong to another group of Culex. To estimate 
N e , we used the method of Xu and Fu 2004 (31) designed 
for highly variable markers such as microsatellites, where 
8=4 N e \i, and u. is the mutation rate. If the parameter 
values needed to get the observed _F ST -values are outside 
what is realistic in this system, we can reject the hypothesis 
of no gene flow between them. 

Results 

The number of alleles, allelic richness, and expected 
heterozygosity did not differ much between the popula- 
tions (Table 1). The global F ST was 0.090 (95% CI: 0.049- 
0.14) using the uncorrected data set and 0.054 (95% 
CI: 0.032-0.80). The pairwise F ST -values were low but 
were not different from zero, except in one case; between 
Denmark and Serbia (Table 2). This result was consistent 
regardless of which data set was used. However, the 
comparisons that involved the samples from France, 
Rumania, and Turkey were highly differentiated from 
the rest of the populations. This can be illustrated in the 
NJ-trees made (Fig. 2). As can be seen in Fig. 2, the 
distance between populations was on average lower using 
the corrected data set. If we exclude the three deviating 
populations, global F ST was now only 0.029 (0.018-0.041) 
and 0.022 (0.015-0.028), respectively, which corresponds 
to G'sT-values of 0.20 and 0.15, respectively. The two 
estimates of F$t was strongly correlated (r=0.99, P< 
0.001; Mantel test). There was a significant correlation 
between geographic and genetic distances [log (distance) 



versus F ST /(1 -F ST ), r=0.31, P =0.0.032 for the cor- 
rected data set, Mantel test]. 

Two of the markers used here were the same that were 
used in Fonseca et al. (4), where two distinct lineages were 
found. In one of the loci CxpGT9, the allele frequency 
spectrum of the three outlying populations (France, 
Romania, and Turkey) was similar to the underground 
populations discovered in Ref. (4) and the allele fre- 
quency spectrum for the rest was close to the one found 
in the above-ground populations found in Ref. (4). The 
difference between the two populations can be seen in 
Fig. 3a, where the allele size 116 [corresponding to allele 
size 114 in Ref. (4)] is shown to be very frequent in 
France, Romania, and Turkey. In Ref. (4), the same allele 
reached a frequency of around 0.5, while in the other 
populations, the frequencies were much lower. Allele size 
124, on the other hand, was not found in Romania, 
Turkey, or France, but in all other populations, at about 
the same frequency. It is also worth noting that allele size 
116, being the most frequent one in the putative under- 
ground populations, was found in almost all of the other 
populations, albeit at a much lower frequency. This 
suggests a gene flow between the underground popula- 
tions and the above-ground populations. The same 
pattern can be seen in locus CQ29 where allele size 178 
reaches a high frequency in the above-ground popula- 
tions, but not in the below-ground populations (Fig. 3b), 
while allele size 184 is very common in France, being 
almost the only allele found, and relatively common in 
Romania, but having a low frequency in Turkey. Both 
alleles are found in both groups, but at a much lower 
frequency, again indicating some gene flow between the 
groups. 

The mean frequency of private alleles was 0.026, which 
when corrected for sample size yields an estimate of the 
number of migrants of 5.74. The number of migrants 
between Malta and Sweden (the most distant popula- 
tions) was as high as 2.04.We found that the estimate of 
effective population, 0 =2.85, which assuming a mutation 
rate for microsatellites (u.) of 5 x 10~ 4 , gives an estimate 
of N e «1500, and with a more conservative estimate of 
10 ~ 4 , gives N c x7l00. To obtain the range of observed 
,F ST -values by drift alone, we would expect N c in the range 
of at least 6000, which is possible given the conservative 



Table 1. Number of individuals sampled, number of alleles (SE), allelic richness, AR fSE), and expected heterozygosity (SE) for the 
different populations 



Denmark Serbia France Croatia Alsace Malta Romania Spain Sweden Turkey 



N 30 30 30 30 30 30 30 30 30 30 



No. alleles 12.1(3.8) 12.8(5.3) 7.4(3.6) 13.0(6.2) 12.5(5.3) 12.5(3.9) 9.5(4.6) 14.0(5.0) 7.1(3.2) 11.6(4.5) 
AR 4.5(0.7) 4.7(0.4) 3.4(1.2) 4.6(0.7) 4.6(0.7) 4.2(0.6) 4.0(1.0) 4.9(0.5) 4.9(0.4) 3.5(4.3) 

H exp 0.85(0.11) 0.87(0.08)0.63(0.27) 0.88(0.08) 0.86(0.09) 0.82(0.10)0.77(0.12) 0.85(0.14)0.90(0.07) 0.68(0.24) 
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Table 2. Pairwise distances with (a) uncorrected F ST and (b) corrected values 
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estimate of mutation, but only at short times since 
separation (Fig. 4). Thus, even if the observed F ST -values 
are not unattainable, the conditions under which this 
might happen are quite restricted. 



Discussion 

The study showed two major results: first, three of the 
analyzed populations (France, Romania, and Turkey) 
were very different from the rest of the populations 
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Fig. 2. Neighbor-joining (NJ) trees illustrating the genetic distance between populations (a) without correction for null alleles, (b) with 
correction, (c) using the lower 2.5 interval, and fd) using the upper 2.5 interval. 



Citation: Infection Ecology and Epidemiology 2012, 2: 12001 - http://dx.doi.org/10.3402/iee.v2i0.12001 



5 

(page number not for citation purpose) 



Mare Lohmus et al. 



1,0- 



s= 0,6- 



0,4 



0,2- 



0,0- 



a) 



IZZ1116 



fl 



Swed Den Alsace Serbia Kroatia Malta Spain France Rum Turkey 
Populations 



1,0 



0,8 



>> 0,6 



cr 

£ 0,4 



0,2 



0,0 



b) 



□ 178 

□ 184 



Swed Den Alsace Serbia Kroatia Malta Spain France Rum Turkey 
Populations 

Fig. 3. Frequencies of two alleles in (a) locus CxpGT9 and (b) QT29 in the different populations. 



(Sweden, Denmark, Alsace, Croatia, Serbia, Spain, and 
Malta) and they also greatly differed from each other 
(Table 2, Fig. 2). In fact, they were so different that there 
are reasons to believe that they belong to differentiated 
lineages. Furthermore, the allele frequency distribution in 
these populations match the one that has been earlier 
found in the underground populations (4). Second, most 
pairwise differences were low and showed a shallow 
structure among the populations apart from the three 
deviating ones (Fig. 2). The low ,F ST -values between 
populations suggest that there is an extensive gene flow 
between the populations, and this was confirmed by the 
simulations where the conditions for differentiation by 
drift alone was far too stringent given what is known 
about these populations. A low F$t by itself does not 
mean that there is extensive gene flow (32) since F ST is 
dependent on the level of heterozygosity. However, even if 



we rescale the values to account for within-population 
heterozygosity, the overall level of differentiation is still 
low in the above-ground populations. Thus, we are 
confident that basically all European populations of the 
above-ground Cx. pipiens share genes. The estimate of 
gene flow was high, which together with the simulations, 
suggests that the genetic similarity among populations is 
due to ongoing gene flow rather than a result of a recent 
shared history and expansion. 

An alternative explanation is that the species only 
very recently spread throughout Europe, and thus due to 
incomplete lineage sorting and subsequent local stochas- 
tic effects, the populations have not diverged more, which 
is reflected in the unclear relationships between the 
populations. However, this species were described by 
Linnaeus in 1758 (33), and there are no reasons to believe 
that it was an uncommon species at that time. Even if the 
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Time since separation EZ1 < 0,025 



Fig. 4. Expected F ST values (shaded areas) in relation to time to separation in generations (x-axis) and effective population size (N e ) 
based on simulations in POWSIM. The line indicates the upper limit of F ST -levels observed in the current study, and thus the observed 
values are only possible under a model of pure genetic drift for parameter values above the line. The dotted line indicates the average of 
the estimated 7V e over all populations. 



species had very recently arrived to Northern Europe 
when Linnaeus described it, 250 years is ample time to 
develop a population structure if gene flow is restricted 
or absent. The results from the power analysis confirm 
this; the expected ,F ST -values are in our simulation higher 
than observed here, even though we used very con- 
servative input-values. This strengthens the support for a 
wide-ranging dispersal of mosquito-individuals across 
large distances in Europe, and hence the possibility of a 
spread of vector-borne pathogens. 

Vector competence of mosquitoes describes their ability 
to transmit pathogens after taking an infected blood meal. 
The ability of vectors transmitting a pathogen is known to 
vary between populations of mosquito species. The exact 
relationship between vector competence and environmen- 
tal factors is however not clear and can be influenced 
by many external factors (8). According to our study, an 
extensive amount of gene exchange, and hence exchange 
of individuals, occurs between populations of Cx. pipiens 
in Europe. Therefore, infectious diseases spread by this 
species may not be associated with a certain vector 
genotype, but rather with suitable environmental condi- 
tions. Consequently, we would expect these pathogens to 
disperse northward with favorable climatic parameters. 

The distribution and intensity of mosquito-borne in- 
fections is strongly influenced by the interactions between 
temperature, vectors, hosts, and pathogen genetics. Tem- 
perature often determines both the latitudinal boundary 
and upper elevation limit of pathogen transmission when 
the extrinsic incubation period becomes greater than the 
life span of the vector (8). Temperature is also linked to 



the intensity and transmission of pathogens. In temperate 
regions, seasonality is one of the key components of the 
climate. In many temperate regions, however, the summer 
temperatures can be as high as temperatures in the tropics. 
The difference with tropics lays in the absence of winters. 
Therefore, tropical mosquito-borne diseases can be intro- 
duced and transmitted during the right season in the 
temperate regions, but in most cases, they will be 
eliminated during the winter (1). Some viruses, however, 
such as the Usutu virus in Austria, have been shown to be 
able to maintain in a natural cycle with a minimum of at 
least 10 days/year with a temperature >30°C (10). It is 
also predicted that Usutu virus will become endemic 
in large parts of Central Europe before the end of the 
century and that optimal environmental conditions for 
outbreaks of this virus will occur in about 10 years from 
now (10, 34). 

A concern associated with pathogens carried by 
Cx. pipiens mosquitoes is the distribution of hybrid 
populations between bf. pipiens and molestus. Because 
pipiens is strictly ornithophilic while molestus has a 
more catholic taste and readily bites humans and other 
mammals, these hybrid populations are thought to be the 
most high-capacity vectors in southern Europe and North 
America (35). This gives rise to concern, should those 
hybrids spread to northern Europe or if another hybridi- 
zation event should occur. It is possible that presently the 
populations of Cx. pipiens in Sweden and the rest of 
northern Europe are less likely to be involved in zoonotic 
transmission of, for example, WNF since they are not 
constituted of hybrids. However, in a recent paper, 
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Reusken et al. (35) showed the presence of hybrid 
pipiens xmolestus mosquitoes in the Netherlands. Since 
our results imply an extensive gene flow between popula- 
tions of Cx. pipiens in Europe, hybrid mosquitoes could 
spread fast over the whole Europe. This could potentially 
alter the epidemiology of many arthropod-borne viruses 
in Northern Europe giving rise to epidemics/epizootics 
similar to the WNF outbreak in North America from 
1999 and onward. Consequently, up-to-date surveillance 
systems are called for all over the Europe as, without 
them, detection and response to emerging vector-borne 
diseases can be severely impaired (36). The emergence of 
WNF in 1999 in New York is a good example of the 
consequences that a delayed response can have on the 
outcome of novel pathogen introduction. Studies show 
that impaired mosquito surveillance can make the man- 
agement costs of epidemics several times higher than 
applying a strategy with sustained surveillance and early 
case detection (36). 
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